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ABSTRACT 



We study the evolution of the broad, double-peaked Ha emission-line profile of the 
LINER/Seyfert 1 nucleus of NGC 1097, using 24 spectra obtained over a time span 
of 11 yrs — from 1991 Nov. through 2002 Oct. While in the first 5 yrs the main 
variation was in the relative intensity of the blue and red peaks, in the last years we 
have also observed an increasing separation between the two peaks, at the same time as 
the integrated flux in the broad line has decreased. We propose a scenario in which the 
emission originates in an asymmetric accretion disk around a supermassive black hole, 
whose source of ionization is getting dimmer, causing the region of maximum emission 
to come closer to the center (and thus to regions of higher projected velocity). We use 
the observations to constrain the evolution of the accretion disk emission and to evaluate 
two models: the elliptical disk model previously found to reproduce the observations 
from 1991 to 1996 and a model of a circular disk with a single spiral arm. In both 
models the peak emissivity of the disk drifts inward with time, while the azimuthal 
orientation of the elliptical disk or the spiral pattern varies with time. In the case of 
the spiral arm model, the whole set of data is consistent with a monotonic precession 
of the spiral pattern, which has completed almost two revolutions since 1991. Thus, 
we favor the spiral arm model, which, through the precession period implies a black 
hole mass that is consistent with the observed stellar velocity dispersion. In contrast, 
the elliptical disk model requires a mass which is an order of magnitude lower. Finally, 
we have found tentative evidence of the emergence of an accretion disk wind, which we 
hope to explore further with future observations. 
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— galaxies: nuclei — galaxies: Seyfert — line: profiles 
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1. Introduction 

The barred spiral galaxy NGC 1097 was the first example of a LINER (low-ionization nuclear 
emission-line region; Heckman 1980) to display broad, double-peaked Balmer lines. The double- 
peaked Ha line (FWHM ^ 7500 km s"^) was discovered by Storchi-Bergmann, Baldwin, &; Wilson 
(1993; hereafter SB93) after a careful subtraction of the starlight from the nuclear spectrum. The 
broad line was not easily recognizable in the original spectrum because its contrast against the 
continuum from the underlying stellar population was fairly low. This discovery was followed 
by other examples, all of them in Hubble Space Telescope (HST) spectra taken through narrow 
apertures in which starlight from the host galaxy is largely excluded. These examples include M81 
(Bower et al. 1996), NGC 4203 (Shields et al. 2000), NGC 4450 (Ho et al. 2000), and NGC 4579 
(Barth et al. 2000). 

Comparison of the double-peaked Ha spectrum of NGC 1097 with earlier observations, as well 
as the later evolution of the double-peaked profile, led to the conclusion that the broad Balmer 
lines were associated with a transient event (e.g., SB93, Eracleous et al. 1995; Storchi-Bergmann 
et al. 1995, 1997; hereafter SB95, SB97), presumably an episode of mass accretion onto the nuclear 
supermassive black hole (SMBH). 

The recent findings of a proportionality between the masses of galactic bulges and that of 
the SMBH in galaxies that are close enough to allow the SMBH's detection via stellar kinematics 
(e.g., Ferrarese & Merritt 2000; Gebhardt et al. 2000) supports their presence in the nucleus of 
most present-day galaxies. Such a SMBH will eventually accrete mass, such as via the capture of 
a star passing closer to the black hole than its tidal radius. If the SMBH has a mass smaller than 
~ lO^M©, a solar-mass star will be tidally disrupted before it is accreted, leading to the formation of 
an accretion disk or ring. As argued by SB95 and SB97, such a process is the preferred interpretation 
for the double-peaked lines of NGC 1097 and may also be applicable to some of the other LINERs 
mentioned above. It is thus possible that we have begun to witness such tidal disruption events, 
which have long been predicted by theory (see review by Rees 1988). 

Monitoring the variability of the double-peaked profiles provides useful constraints on the above 
scenario. Therefore, we have been following the evolution of the Ha profile of NGC 1097 since it 
was first discovered in 1991. In this paper we present the result of 11 years of monitoring and use 
the profile observations to constrain and refine scenarios for their origin, such as the precessing 
elliptical ring suggested by Eracleous et al. (1995) and SB95. 

In §2 we describe the observations and analysis of the data, in §3 we present the results, and in 
§4 wc discuss their interpretation and present two possible accretion disk models to reproduce the 
profile variations. We explore in §5 the significance of our findings and compare them with those 
for other objects. The conclusions of this work are presented in §6. 
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2. Observations 

We observed the nuclear spectrum of NGC 1097 from 1991 to 2002, up to several times a year, 
using mostly the CTIO 4m Blanco and 1.5m telescopes, and the KPNO 2.1m telescope. We also 
obtained a few spectra with the Keck II, Lick 3m, MDM 2.4m, and ESO NTT telescopes. The log 
of observations is given in Table 1, where we also list the slit and extraction window widths. 

The first observations (SB93, SB95, SB97) were obtained through a 2" slit and an extraction 
window of 4". We later used a narrower slit of 1'.'5 for the observations with the large telescopes 
(3.6m and 4m) in order to improve the contrast between the double-peaked profile and the starlight 
spectrum, but keeping it wide enough to avoid spurious variations due to seeing differences. We 
have also kept the extraction window large enough to allow the normalization of the double-peaked 
broad line flux using the fluxes of the narrow emission lines, as discussed below. 

In order to isolate the double-peaked profile, we subtracted the contribution of the stellar 
population, as described in SB93, SB95, and SB97. In the case of the observations with the Blanco 
and NTT telescopes, the stellar population spectrum was extracted from windows adjacent to the 
nuclear window, typically from regions between 3" and 6" from the nucleus. The nuclear spectrum 
is redder than that of the stellar population at these windows. Attributing the difference in color 
to reddening by dust, we estimate a value of E{B — V) ^ 0.05 mag. We have thus reddened 
the stellar population spectra accordingly and normalized them to the flux of the continuum in 
the nuclear spectrum. In addition, we have excized the residual narrow-line emission from the 
stellar population spectra (using as a reference spectra of disk galaxies free of emission lines), and 
have subtracted it from the nuclear spectrum. The spectra obtained with the smaller telescopes 
(CTIO 1.5m and KPNO 2.1m) have a signal-to-noise ratio too low outside the nucleus to allow 
the extraction of a population template; for these we used the starlight spectra extracted from the 
CTIO 4m observations to subtract the stellar population contribution. The same was done for the 
Keck and MDM spectra. The process of starlight subtraction is more important in the most recent 
spectra because the double-peaked line is fading (see discussion below) and its contrast to the 
stellar population spcctrTmi is decreasing. Thus a careful match of the stellar population spectrum 
is necessary before the subtraction. 

The last step was to normalize the broad-line flux, which we have done using the fluxes of 
the narrow emission lines Ha, [N II] AA6548, 6584, and [S II] AA6717, 6731, under the assumption 
that they have not varied during the 11 years of observations. This is a reasonable assumption, 
considering that the average extraction windows correspond to ~ 150 x 300 pc^ at the galaxy. 
For the normalization we used as reference the spectrum of 1996 Jan. 24, as follows. For each 
spectrum we determined the average ratio of the [N II] -|- Ha(narrow) and [S II] fluxes with respect 
to the 1996 Jan. 24 spectrum and multiplied the spectrum by this average ratio. This scaling 
factor is listed for each spectrum in the last column of Table 1. The success of this procedure is 
illustrated in Figure 1, which shows that profiles obtained approximately at the same epoch (within 
two months of each other) using different telescopes match each other very well after the above 
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normalization. This result also shows that the small differences in apertures of the different spectra 
does not significantly affect the normalization. 

The Keck spectra are the only exceptions, as we had to "re-normalize" them to ~80% of the 
resulting flux after the normalization described above, to allow an agreement with the spectra 
obtained at approximately the same epoch with the other large telescopes. We attribute this result 
to the narrower slit used in the Keck spectra (!")) suggesting that the narrow- line emission extends 
beyond 1" although not much more, as the normalization works for the larger apertures. The 
scaling factors listed in the last column of Table 1 for the Keck spectra are already corrected for 
this additional normalization. 

After normalization, Figure 1 shows that there is reasonable agreement between the CTIO 4m 
and KPNO 2.1m spectra of 1997 Sep. There is excellent agreement between the spectra with the 
best signal-to-noise ratios: Keck (thin continuous line; epochs 1998 Jan. and 1999 Jan.) and CTIO 
4m (thick continuous lines; epochs 1998 Feb. and 1998 Dec). In comparison, the MDM spectrum 
of 1998 Dec. shows a broad dip around 6530 A which is not present in the Keck and CTIO spectra. 
There is also good agreement between the Keck spectra of 2000 Oct. and the KPNO 2.1m spectra 
of 2000 Sep. The spectra obtained with the CTIO 1.5m have the poorest signal-to-noise ratios: the 
profile from 1998 Jan. has very little flux and is not compatible with the profiles measured both 
with the Keck and CTIO 4m at approximately the same epoch. On the other hand the 1.5m profile 
of 1998 Oct. agrees with that of the KPNO one of 1998 Sep. and with the CTIO 4m one of 1998 
Dec. 

On the basis of Figure 1 we conclude that we should give more weight in our analysis to the 
results based on the Keck, CTIO 4m, and NTT 3.6m spectra relative to the ones obtained with 
the smaller telescopes. This figure also illustrates that changes in the shape of the double-peaked 
profiles seem to occur on time scales longer than a few (^ 3-4) months. 

3. Results 

In Figures 2 and 3 we show a sequence of selected Ha profiles, from 1991 Nov. 2 to 2002 Oct. 
10. Our criteria for the selection was to obtain the most complete time coverage using the available 
data, but using spectra having the highest signal-to-noise ratios in cases of multiple spectra at a 
given epoch, as discussed above. 

During the evolution of the double-peaked Ha profile shown in Figures 2 and 3, the most 
noticeable changes occur in the relative strengths of the blue and red peaks. The profile starts 
with the red peak stronger than the blue in 1991 Nov. and evolves to similar strengths of the two 
peaks in 1994 Jan., then to the blue peak stronger than the red in 1996 Jan. The spectrum of 1997 
Sep. shows again the red peak stronger than the blue, then the two peaks show similar strengths 
by mid 1998, and then the blue peak is stronger than the red by 1999. In the observation of 2001 
Sep. the blue and red peaks have similar strengths, and in 2002 Oct. the blue peak is stronger 
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again. Another important change that could only be noticed after accumulating data over so many- 
years is a decrease of the broad-line flux at the same time as the blue and red peaks are moving 
progressively farther apart, as indicated by the measurements below. 

In order to quantify the observed variations, we measured the wavelengths of the blue and 
red peaks, their peak flux densities, and the flux of the broad line, which we list in Tabic 2. The 
wavelengths of the blue and red peaks {Xb and Xr) and corresponding peak flux densities (Fg and 
Fr) were obtained by fitting Gaussians separately to the blue and red peaks of the profiles and 
measuring the central wavelengths and peak fluxes. These measurements were only possible in the 
large-telescope spectra, as the fits to the blue and red peaks did not consistently converge to a imique 
central wavelength in spectra with low signal-to-noise ratio. The flux of the broad line (-Fferoad) was 
measured by integrating the flux of the broad plus the narrow emission lines and then subtracting 
the fluxes of the narrow lines. Uncertainties in the wavelengths and fluxes were estimated by 
repeating the fits several times after changing by a few angstroms the limits of the windows used 
to fit the Gaussians and to measure the broad-line flux. The resulting average uncertainties are 
~ 5 A for Xb and Xr, 5 x 10~^^ erg cm" "2 s-^ A-^ for Fb and Fr, and ~ 2 x lO"^^ erg cm s 

for F^yroad- 

Figure 4 shows the temporal variations of the blue and red peak velocities {Vb and Vr), 
obtained from A^, Xr, and the average wavelength of the narrow Ha line as reference (6591.7 A), 
as well as the temporal variation of Vr — Vb and of Ff^road- We observe a blueshift of the blue peak 
from ~ —2800 km in the first years to ~ —4500 km s~^ in 2001, and a redshift of the red peak 
from ~ 3500 km to ~ 4500 km s~^ over the same period. The resulting Vr — Vr increases from 
^ 6000 km s"^ to ~ 9000 km s~^. In Figure 5 we show the temporal variation of Fr/Fr. 

4. Interpretation and Comparison With Models 

4.1. Preliminary Inferences from the Data 

The continued variation of Fr / Fr and, in particular, the strong red peak observed in 1997 Sep. 
still support the origin of the profile in an accretion disk that is not axisymmetric. The asymmetry 
could be the result of an eccentricity (e.g., Eracleous et al. 1995; SB95), or a circular disk with a 
disturbance, such as a bright spot (Newman et al. 1997; Zheng, Veilleux, & Grandi 1991) or spiral 
arm (e.g., Chakrabarti &: Wiita 1993; Gilbert ct al. 1999). The precession of the elliptical disk or 
of the disturbance in the circular disk would be the origin of the variation of Fr / Fr with time. 

The increasing separation of the blue and red peaks (§3) indicates that the bulk of the broad-line 
emission is coming from regions of successively larger velocities. In the accretion disk scenario, this 
means that the bulk of emission is coming from successively smaller radii. The above result could 
suggest at first glance that we are witnessing the infall of matter toward the SMBH. Unfortunately, 
this seems not to be possible over ten years of observations. A noticeable decrease in radius should 
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occur in the viscous time scale, which should be larger than 10'^ years in the case of NGC 1097 
(Prank, King, k, Raine 1992, p. 99, using the black hole mass obtained in Sec. 5.2). 

An alternative possibility is borne out by the observation that, at the same time that the profile 
is getting broader, the line flux is decreasing (in 2001 Sep. the broad Ha flux is ~l/3 of the flux 
of 1991 Nov.), which suggests that the ionizing source is getting dimmer. Thus, a more plausible 
interpretation is that the dimming of the ionizing source is causing the region of maximum emission 
to come closer to the center, to regions of larger velocities and thus giving rise to broader proflles. 

Inspection of Figure 4 shows that there seems indeed to be an inverse correlation between 
Vr — Vb and Fbroad, not only as a general trend along the 11 years covered by the observations, 
but also from epoch to epoch. For example, in 1996 Jan. the flux of the line showed an increase 
relative to the two previous epochs, while Vr — Vb decreased. The same seems to have happened in 
Sep.-Oct. 2000. This inverse correlation can be better observed in Figure 6, where we plot Vr — Vb 
against Fhroad- 

The above effect (the inverse correlation between the ionizing flux and width of the profile) has 
been previously discussed by Dumont &: Collin-Souffrin (1990). These authors explain the existence 
of a "radius for maximum emission" or "saturation radius" as follows. In the line-emitting region of 
the accretion disk there are two effects contributing to the line intensity: the increase of the emitting 
surface with distance from the center, and the decrease of the intensity of ionizing radiation with 
radius. The ionizing flux can become quite high close to the center, leading to ionization of the 
disk and saturation of the line emission (cooling is mainly due to free-bound and free-free processes 
and not to lines). The radius at which this saturation sets in is the radius of maximum emission. 

Finally, we notice that in Figure 4, the velocity of the blue peak shows a steeper variation 
(blueshift) than that of the red peak (redshift), implying a net blueshift of the center between 
the two peaks of ~ —450 km between 1991 Nov. and 2001 Sep. One possibility is that this 
blueshift is due to a wind originating in the disk, which would be consistent with observations of 
Fe II absorption lines in a recent HST spectrum of NGC 1097 (Eracleous 2002) which indicate the 
presence of an outflow. 

With these considerations in mind, we have tried to fit the variable line profiles with two 
different models, as we describe in detail below. Both models rely on a non-axisymmetric disk 
structure and its precession to reproduce the changes in the relative strengths of the two peaks. 
However, this mechanism cannot explain the increasing separation of the two peaks. The latter 
effect requires us to change the radius of maximum line emission in the disk by adjusting the 
prescription for the disk emissivity. 
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4.2. Elliptical Disk Model 

In Eracleous et al. (1995) we showed that an elhptical accretion ring model could reproduce 
the broad double-peaked Ha profile of NGC 1097, as observed in 1991 Nov. After following the 
variations of the line profile from 1991 to 1996 (SB97), we developed the model further so that the 
plane of the ring was inclined at i = 34° to the line of sight, and its line-emitting part corresponded 
to a range of pericenter distances between = 1300 and ^2 = 1600, where ^ is expressed in units 
of the gravitational radius rg = GMjS . The ring was circular up to a radius = 1400, and from 
there the eccentricity increased linearly from e = to 0.45. The ring had a broadening parameter 
G = 1200 km s~^ and an emissivity varying as e(^) = eo The variation of the observed line 
profiles from 1991 to 1996 were reproduced by the precession of the ring with a period of of 25 
years, which led to an estimate of the central black hole mass of 9 x 10^ M0. 

We tried to use the same elliptical ring model to reproduce all the double-peaked Ha profiles 
observed from 1991 to 2002, by changing only the orientation angle, 4>o, under the assumption that 
the precession continued. However, we were not able to reproduce the line profiles observed after 
1996; continued precession alone could explain the variation in the relative heights of the two peaks, 
but not the observed velocity shifts of the blue and red peaks. After exploring the effect of other 
parameters on the shape of the line profile, we found that only changes in the inner pericenter 
distance of the ring (from ^1 = 1300 to 450) could lead to the observed increase in the separation of 
the two peaks, Xr — Xb- Such a dramatic change in the physical location of a given volume of gas 
could not occur on a time scale of 5 years, however, as we noted in the previous section. Rather, 
the emission properties of the ring must have changed so that the radius of maximum emission 
became smaller. Therefore, to produce a self-consistent and physically plausible scenario for the 
entire set of observations, we were led to revise the model parameters. This revision represents a 
refinement of earlier models, which is possible only after following the variability of the line profile 
over a baseline of 10 years. We now propose that the ring has always had a smaller inner radius of 
^1 = 450 (and now resembles a disk more than a ring), and that the increasing width of the profile 
is due to a change in emissivity such that, in more recent epochs, the inner regions of the disk make 
a more significant relative contribution to the total line emission. This picture is also consistent 
with the observational result that the flux of the broad line is decreasing, as discussed above. 

To explain the increasing separation between the two peaks, we experimented with several 
emissivity laws. We found that we could best reproduce the observations with a broken power 
law, such that inside a radius the power-law index is gi = — 1 and outside this radius q2 = 1. 
This means that the emissivity increases with radius from ^1 to and then decreases with radius 
beyond ^q. These q values are also physically motivated: in the "saturation" region, the increase 
in emissivity is due to the increase in the area of the ring, while beyond ^g, the decrease can be 
understood as due to the radiation field dropping off as r^^ and the increase of the emitting area 
as r^. 

In this parameterization of the disk emissivity, varying the "break radius" allows us to 
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reproduce the shifts of the blue and red peaks. In particular, allowing the "break radius" to 

decrease with time shifts the region of maximum emission closer to the center and results in an 
increase in the separation between the two peaks. At the same time, since the relative heights of 
the two peaks have continued to change, we allow the disk to precess (as in SB97). 

To fit individual line profiles we fixed the inner and outer pericenter distances of the line- 
emitting disk to ^1 = 450 and ^2 = 1600 (respectively), the inclination angle to i = 34°, and the 
broadening parameter to cr = 1200 km s~^. The eccentricity was allowed to increase linearly with 
radius, from e(^i) = to 6(^2) = 0.45. The emissivity of the disk (e oc was set to increase with 
radius in the inner disk (i.e., q'l = — 1 up to a radius ^q) and then decrease (^2 = 1)- The orientation 
angle of the disk, (po, and the emissivity break radius, $,q, were the only parameters allowed to vary 
to represent a precession of the disk and a decline in the illumination. We first searched for models 
that could reproduce the primary observables, namely the velocities and relative strengths of the 
two peaks. The best-fitting models resulting from the above procedure are compared with the 
observations in Figures 7 and 8, and the model parameters are summarized in Table 3. 

In the last epochs, beginning in 1998 Dec, the fits of the model to the data show some 
improvement if we allow a small blueshift of the central wavelength of the profile (~10A) relative 
to the wavelength of the narrow Ha line (6591.7 A — which we have so far adopted as the center 
of the broad line in the models). We show these improved fits as dotted lines in Figure 8. This 
blueshift of the broad line was already detected in the observed evolution of the blue and red peak 
velocities (Fig. 4), and was discussed at the end of the previous subsection. 

The evolution of and (pQ arc shown in Figure 9. There is very little change in until 
1996, but then in 1997 it begins to decrease, causing the relative contribution of the inner parts of 
the disk to become more important and producing the broadening of the profile. The orientation 
angle of the disk, represented by 00, shows a small range of values, 131 < 00 < 209, as listed in 
Table 3. Noticing that changes in orientation by half a revolution result in very similar (although 
not identical) model profiles, we have exploited this property in order to construct a series of models 
that would approximately represent a precession of the disk. The best compromise we could obtain 
between a good fit of the profiles and the precession scenario is illustrated in the bottom panel 
of Figure 9. Assuming that the precession rate of the disk is constant, we carry out a linear 
least-squares fit to the precession phase shown in Figure 9. We adopt a functional form given by 

<f>0 _ t - to /-.X 
^ prec 

where t is the time in years, to is an arbitrary zero point of the precession cycle, and Pprec is 
the precession period. The best fit, shown as a solid line in Figure 9, yields to = 1990.8 and 
^prcc = 7.1 years, which implies that the disk has completed almost two revolutions since the epoch 
of the first observation as one would expect from the behavior of the peak strength ratio, Fb/Fr, 
shown in Figure 5. In Figure 10 we show a cartoon of the accretion disk in two epochs, the first 
(1991.8) and last (2001.7) fitted by the model, illustrating the change in the disk emissivity in 10 
years. 
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in conclusion, we note that even though models that represent a smooth precession of the 
eccentric disk are preferable on physical grounds, they do not provide as good a description of the 
line profiles as the models whose parameters are listed in Table 3. 



4.3. Circular Disk Plus Spiral Arm Model 

An alternative to the elliptical disk model is a model that invokes a non-axisymmetric per- 
turbation of the emissivity of a circular disk. Models of this type were applied to the variable Ha 
profile of 3C 390.3 by Zheng, Veilleux, & Grandi (1991; a "bright sector" in the disk) and by Gilbert 
et al. (1999; a single-arm spiral, inspired by the work of Chakrabarti & Wiita 1993, 1994). An 
important feature of these models is that the perturbation spans a relatively wide range in radius 
(and projected velocity), which leads to a modulation in the relative strengths of the two peaks of 
the line profile without creating a "third peak" that moves across the profile. Spiral arm models 
have additional appeal because spiral arms are known to exist in other astrophysical disks, such 
as the disks of spiral galaxies and the accretion disks of some cataclysmic variables (e.g., Steeghs, 
Harlaftis, k, Horne 1997; Baptista & Catalan 2000; Patterson, Halpern, & Shambrook 1993). 

We adopt the formulation of Gilbert et al. (1999), according to which the line-emitting portion 
of the disk is circular, has an inclination angle z, and is bounded by inner and outer radii and ^2- 
Superposed on the axisymmetric emissivity pattern (e oc ^~^) there is an emissivity perturbation 
which has the form of a single spiral arm. We adopted a single-arm pattern after experimenting 
with patterns having two or more spiral arms. We found that two or more spiral arms could not 
reproduce the changing asymmetry of the Ha line profile of NGC 1097, a conclusion also reached 
by Gilbert et al. (1999) for 3C 390.3. 



The prescription for the total emissivity of the disk is given by 



e(C,</') = eor^ <{ 1 + ^exp 



A 

+ -exp 



.__(27r-</, + V'o)' 



(2) 



where the two exponentials represent the decay of the emissivity of the spiral arm with azimuthal 
distance on either side of the ridge line. The parameter 6 represents the azimuthal width of the 
spiral arm (FWHM), while the quantity (f> — tpo is the azimuthal distance from the ridge of the 
spiral arm, which is defined by 

Vo = (Po-\ — , (3) 

tanp 

where p is the pitch angle and ^sp is the innermost radius of the spiral arm (also used as a fidu- 
cial radius for normalization; the spiral arm is set to begin at the outer radius of the disk, ^2)- 
The parameter A represents the brightness contrast between the spiral arm and the underlying, 
axisymmetric disk, while sets the azimuthal orientation of the spiral pattern. 



To fit individual line profiles, we followed the procedure described in §4.2. We found that we 
could keep the inner and outer radii of the disk at the values of the inner and outer pericenter 
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distances of the elliptical disk model, namely = 450 and ^2 = 1600. We also kept the inclination 
at i = 34° and the broadening parameter at o" = 1200 km s~^. The only parameters we needed to 
vary in order to reproduce the evolution of the profiles were q, A, and (pQ. The remaining spiral 
arm parameters were kept at the following values: ^sp = ^i, P = —50°, and 6 = 70°. As in the case 
of the elliptical disk model, the variation of cpQ can explain the modulation of the relative strengths 
of the two peaks of the Ha line but not the increasing separation between them. Therefore, in the 
spirit of the discussion of §4.1, we allowed the emissivity power-law index to vary so as to shift the 
region of maximum emission to smaller radii. 

The best-fitting models are compared to the observed line profiles in Figures 11 and 12, while 
the corresponding parameter values (for A, q, and ^o) are listed in Table 4. The spiral arm contrast, 
parameterized by A, has decreased from 3 to 2 since the first observations, while q has increased, 
favoring the emission of the inner part of the disk relative to that of the outer parts in the later 
epochs and producing the broadening of the profile. The 4>o values were constrained to allow a 
monotonic precession of the spiral pattern, as illustrated in Figure 13, where we show how q and 
(f)o evolve with time. The latter figure shows that the spiral model accomodates better a monotonic 
precession (of the spiral pattern) than the elliptical model. The solid line in the lower panel of Figure 
13 is a fit to the precession phase following Eq. (1), yielding to = 1991.3 and Pprcc = 5.5±0.3 years. 
Such a period implies, again, that the disk has completed almost two revolutions since the epoch 
of the first observation as one would expect from the behavior of the peak strength ratio {Fb/Fr) 
shown in Figure 5. 

In Figure 14 we illustrate a cartoon of the disk in the first and last epochs fitted by the model. 
Notice that the spiral arm is quite "broad" and hard to observe in the last epoch due to the low 
contrast relative to the disk as implied by the best-fitting parameters (Table 4). 

5. Discussion 

5.1. Comparison With Studies of Other Objects 

The variation of the relative strengths of the two peaks of the Ha line of NGC 1097 appears to 
follow the same trend as in broad-line radio galaxies. Examples include 3C 390.3 (Zheng, Veilleux, 
& Grandi 1991; Gilbert et al. 1999), 3C 332 (Gilbert et al. 1999), Pictor A (Eracleous & Halpern 
1998), and Arp 102B (Newman et al. 1997). On the other hand, drifts in the peak wavelengths 
of double-peaked profiles are seldom observed. One case previously reported is that of the radio 
galaxy 3C 390.3, for which Gaskell (1996), Eracleous et al. (1997), and Gilbert et al. (1999) found 
a drift in the velocity of the blue peak between 1970 and 1996. But in 3C 390.3, Eracleous et al. 
(1997) found that the two peaks moved together, unlike what is observed in the case of NGC 1097, 
in which the two peaks move apart from each other. 

On the other hand, Shapovalova et al. (2001), using observations of the H/3 profile of 3C 390.3 
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from 1995 to 2000 (later than those reported by Eracleous et al. 1997), claim to have found drifts 
of the two peaks in opposite directions (or velocities), arguing that the difference between the 
velocities of the blue and red peaks is inversely correlated with the continuum flux. They note 
that when the H/3 flux decreases, the radial velocity difference between the red and blue peaks 
increase, while when the H/3 flux increases, the velocity difference between the two peaks decreases. 
We point out, however, that the measurement of the wavelength of the red peak of the H/? line in 
3C 390.3 is quite uncertain because this peak is severely contaminated by the [O III] AA 4959, 5007 
emission lines. 

The behavior of the 3C 390.3 broad H/3 profile as described by Shapovalova et al. (2001) would 
be similar to the case of NGC 1097 reported here (as discussed in Sec.4.1), and they have proposed 
a similar interpretation: the flux increase of H/3 is due to the increase in the source luminosity, 
which thus ionizes regions farther from the nucleus (lower velocity); the flux decrease of H/3 is due 
to a decrease in the ionizing luminosity, which shifts the region of maximum line emission closer to 
the center of the disk (higher velocity) . 



5.2. Implications of the Observational Results and Models 

The pattern precession period inferred from the model fits can be connected to the mass of 
the central black hole. Although eccentric disks and spiral arms are related structures (see, for 
example, the illustration in Adams, Ruden, h Shu 1989), in our application the precession rates 
are set by different mechanisms. 

In the case of an eccentric disk, the precession rate is controlled by the relativistic advance of 

the pericenter. Therefore, we adapt the formula given by Weinberg (1972, p. 197), which we write 

as ^ 

27r 1 + e GM .in 1 + e ,./2 

^Prec = y (T^^ ^ ^ ' = (13^ ' 

where Mq is the mass of the black hole in units of 10^ and ^3 is the pericenter distance of the 
orbit in units of 10^ rg. If the outer, most eccentric ring of the disk is precessing as prescribed by 
the advance of pericenter, then ^3 = 1.6, e = 0.45, and, for a period of 7.1 years obtained in §4.2, 
the resulting black hole mass is 6 x 10^ M©. If the inner rings impose a faster precession rate, the 
resulting mass will be somewhat larger, with an upper limit corresponding to the innermost radius 
of ^3 = 0-45 and e = 0, which yield a black hole mass of 5 x 10^ Mq. 

In the case of the precessing spiral pattern, we can connect the precession period to the black 
hole mass using the results of Laughlin & Korchagin (1996), who find that the pattern period is 

several times to an order of magnitude longer than the dynamical time. An extreme upper bound 
to the pattern period is the sound-crossing time of the disk. The dynamical and sound-crossing 
time scales are respectively given by 

Tdyn = 2 Me days (5) 
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and Ts = 8 Me ^ T5 ' months (6) 

(adapted from Prank et al. 1992), where now ^3 is the radius in the disk in units of 10^ rg and 
is the temperature in units of 10^ K. By comparing these time scales with the precession period of 
5.5 years inferred in §4.3, we estimate a black hole mass in the range 10^-10^ M©. 

Another scenario discussed briefly in SB97 involves the precession of a radiation-induced warp 
in the disk (Pringle 1996). In Eq. (4) of SB97 we estimated a precession period of the order of a 
decade assuming a black hole mass of 10^ Mq and a mass of the accretion disk of 10^"^ Mq. Since 
this disk mass is, in fact, a lower limit (see SB97), the precession period determined earlier implies 
an upper limit on the black hole mass of M < 10^ Mq in the context of this picture. 

All of the above scenarios can be assessed using a black hole mass derived from the velocity 
dispersion of the stars in the nucleus of NGC 1097. Emsellem et al. (2001) measure a velocity 
dispersion of 160±20 km s~^ in the inner 300 pc of the bulge of NGC 1097. Using the mass- velocity 
dispersion relation of Tremaine et al. (2002), we infer a black hole mass in the range (2-9) X 10^ M0, 
where the uncertainty is dominated by the dispersion about the best-fitting relation (see Figure 
8 of Tremaine et al.). This black hole mass disfavors both the elliptical disk and precessing warp 
scenarios, leaving the spiral arm scenario as the best candidate. The spiral arm scenario is also 
favored by the fact that it provides the best description of the variable line profiles, by better 
"accommodating" the precession scenario, as described in §4. 



6. Summary, Conclusions, and Open Questions 

We have used spectroscopic observations of the broad double-peaked Ho; profile of the nucleus 
of NGC 1097 spanning the period between 1991 Nov. and 2002 Oct. in order to constrain and refine 
models for the variability of the line profile. The profiles display the following types of variations 
over these 11 years. 

1. A smooth, long-term variation in the relative strengths of the blue and red peaks. 

2. A decrease in wavelength of the blue peak and a corresponding increase in wavelength of the 
red peak — in other words, a broadening of the double-peaked profile, suggesting that the 
bulk of the emission is shifting inward to gas that is moving at higher velocities. 

3. A decrease in the integrated luminosity of the line, so that in 2001 Sep. it was approximately 
1/3 of its value of 1991 Nov. 

4. An inverse correlation between the width and luminosity of the line. 

5. A shift of the central wavelength of the profile (defined as the average wavelength of the 
two peaks) from approximately zero velocity relative to the narrow lines in 1991 to about 
-450 km s"^ in 2001. 
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The evolution of the double-peaked Ha profile of NGC 1097, described above, lends additional 

support to our earlier suggestion that the variations are a result of a precessing, non-axisymmetric 
pattern in an accretion disk (Eracleous et al. 1995; SB95; SB97). We have compared the observed 
line profiles with models of an eccentric (elliptical) disk and a disk with a spiral arm and have 
shown that both models can reproduce the data. However, we favor the spiral arm model because 
(1) it can reproduce the evolution of the double-peaked profile as due to the precession of the spiral 
arm, while the elliptical disk model cannot reproduce a smooth precession; and (2) it implies a 
black hole mass that is consistent with what is inferred from the radial velocity dispersion of the 
stars in the nucleus, while the elliptical disk model requires a mass which is an order of magnitude 
lower. 

We have further found that the dimming of the broad double-peaked Ha line is coupled with 
the increasing separation of the blue and red peaks. We have interpreted this behavior as a result of 
the dimming of the source of ionizing radiation, which makes the maximum emission region of the 
disk to drift inward, and thus to regions of higher orbital velocities. This dimming of the central 
source is progressing independently of dynamical phenomena in the disk, which is a reasonable 
outcome of the tidal disruption scenario that explains the abrupt appearance of the double-peaked 
lines. The coupling between the luminosity and width of the line seems to hold also at shorter time 
scales: an increase in the flux of the line is associated with a decrease in the line width, while a 
decrease in flux is associated with an increase in line width. 

In the process of studying the evolution of the Ha profile and fitting disk models to them we 
found a gradual blueshift, especially at later epochs. This blueshift manifests itself as a gradual 
shift of the average wavelength of the two peaks. We have not attempted to model this effect since 
our data so far cannot yield strong constraints. Nevertheless, we speculate that this blueshift may 
be associated with the onset of an accretion disk wind (Elvis 2000). Outflows with typical velocities 
of 1000 km s^^ — consistent with the blueshift of '--^400-500 km s^^ we have measured for NGC 1097 
if we consider the inclination of 34° for the disk — have been frequently observed in the spectra of 
Seyfert galaxies as absorption lines in the UV (e.g., Kraemer, Crenshaw, & Gabel 2001) and X-rays 
(the "warm absorber"), as well as in the emission lines of the narrow- line region. Models for these 
outflows attribute their origin in the accretion disk, as, for example, a hydromagnetically driven 
wind of discrete clouds (Emmering, Blandford, &: Shlosman 1992), or radiatively driven continuous 
winds (Murray et al. 1995; Proga, Stone, &: Kallman 2000). The presence of an outflow in the 
case of NGC 1097 is also consistent with observations of Fe 11 absorption lines in a recent HST UV 
spectrum (Eracleous 2002). Future observations of the variability of the Balmer-line profiles should 
shed more light on the issue of the possible wind. 
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Table 1. Observation Log 



Date 




Exposure 


Aperture 


Scaling 


(mm/dd/yy) 


Telescope 


(s) 


(") 


factor 


11/02/91 


CTIO 4m 


2700 


2.0x4.0 


1.12 


10/05/92 


CTIO 4m 


685 


2.0x4.0 


1.32 


09/12/93 


Lick 3m 


1200 


2.0x4.0 


0.52 


01/05/94 


CTIO 4m 


1800 


2.0x4.0 


1.14 


12/06/94 


CTIO 4m 


1200 


1.5x3.7 


1.46 


01/24/96 


CTIO 4m 


1200 


2.0x4.0 


1.00 


10/10/96 


Lick 3m 


3600 


2.0x4.7 


0.91 


01/02/97 


CTIO 1.5m 


3600 


1.8x2.6 


0.66 


09/24/97 


CTIO 4m 


1800 


1.5x3.7 


1.49 


09/27/97 


KPNO 2.1m 


3600 


1.9x3.8 


1.22 


01/02/98 


CTIO 1.5m 


6000 


1.8x2.6 


1.79 


01/17/98 


Keck II 


500 


1.0x4.8 


1.28 


02/01/98 


CTIO 4m 


1800 


1.5x3.7 


1.27 


10/15/98 


KPNO 2.1m 


5400 


1.9x4.6 


1.42 


10/20/98 


CTIO 1.5m 


5400 


1.8x6.5 


0.87 


12/16/98 


MDM 2.4m 


2400 


1.5x2.9 


1.19 


12/28/98 


CTIO 4m 


1800 


1.5x3.7 


1.29 


01/06/99 


Keck II 


300 


1.0x2.4 


1.11 


12/02/99 


KPNO 2.1m 


3600 


1.9x3.0 


1.21 


09/23/00 


KPNO 2.1m 


3600 


1.9x3.0 


1.10 


10/11/00 


Keck II 


300 


1.0x2.4 


0.70 


03/01/01 


CTIO 4m 


900 


1.5x3.7 


1.22 


09/14/01 


NTT 3.6m 


1800 


1.5x3.7 


1.33 


10/10/02 


KPNO 2.1m 


3600 


1.9x3.0 


1.17 
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Table 2. Measured Properties 



Epoch 


Ab(A) 






Ai?(A) 














Fr^ 




Fbroad 


1991. 


,84 


6528.0± 


3 


.0 


6670.4± 


1. 


,9 


0, 


,754± 


0, 


,020 


1. 


,172± 


0, 


.020 


227.0± 9.6 


1992. 


.76 


6535.9± 


2, 


.4 


6664.6± 


2. 


,7 


0, 


.801± 


0, 


.030 


1. 


,122± 


0, 


.003 


204.3± 13.1 


1994. 


,01 


6530.9± 


0, 


.4 


6671. 1± 


2. 


,2 


0, 


.497± 


0, 


.010 


0. 


,564± 


0, 


.036 


120.3± 6.1 


1994. 


,92 


6519.3± 


2, 


.2 


6676.8± 


2. 


,2 


0, 


.429± 


0, 


.024 


0. 


,368± 


0, 


.032 


lOl.Oib 14.0 


1996, 


,06 


6528. 1± 





.5 


6672.3± 


3. 


,3 


0, 


,854± 


0, 


,002 


0. 


,628± 





.020 


165. 7± 3.1 


1997. 


,73 


6508.6± 


5, 


.5 


6675.7± 


2. 


,7 


0, 


,304± 


0, 


,017 


0. 


,545± 


0, 


.016 


102. 1± 3.8 


1998. 


,07 


6504.7± 


4, 


.4 


6679.3± 


1. 


,0 


0, 


,398± 


0, 


,005 


0. 


,430± 


0, 


.006 


102.3± 11.1 


1998. 


,09 


6510.1± 


2, 


.5 


6679.4± 


0. 


,7 


0, 


.335± 


0, 


.017 


0. 


,461± 


0, 


.007 


104.8± 1.8 


1998. 


,80 


6502.1± 


2, 


.5 


6686.5± 


6. 


,5 


0, 


.334± 


0, 


.015 


0. 


,292± 


0, 


.012 


76.6± 6.0 


1998. 


,99 


6496.6± 


2, 


.3 


6690.5± 


5. 


,8 


0, 


.283± 


0, 


.009 


0. 


,225± 


0, 


.002 


70.1± 2.1 


1999. 


,02 


6494.5± 


2, 


.2 


6687.2± 


7. 


,5 


0, 


,302± 


0, 


,010 


0. 


,227± 





.006 


71.8± 12.0 


2000. 


,73 


6491. 8± 


2 


.5 


6673.5± 


4. 


,0 


0, 


,502± 


0, 


,020 


0. 


,380± 


0, 


.015 


102.3± 11.0 


2000. 


,79 


6493.1± 


2, 


.3 


6674.4± 


3. 


,2 


0, 


,441± 


0, 


,021 


0. 


,298± 


0, 


.013 


lOO.Oib 7.5 


2001. 


,16 


6495.1± 


1, 


.6 


6687.9± 


2. 


,8 


0, 


.338± 


0, 


.001 


0. 


,254± 


0, 


.005 


80.2± 6.8 


2001. 


,70 


6496.4± 


2, 


.0 


6691.4± 


1. 


,1 


0, 


.326± 


0, 


.007 


0. 


,338± 


0, 


.012 


84.4± 8.5 



^Peak flux density (10 erg cm ^ s ^ A 
^Integrated flux (10~^^ erg cm~^ 
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Table 3. Elliptical Disk Model Parameters 







6n(°) 


1991.84 


1260 


152 


1992.76 


1260 


152 


1994.01 


1200 


171 


1994.92 


1310 


182 


1996.06 


1260 


189 


1997.73 


1140 


131 


1998.09 


1110 


164 


1998.99 


940 


202 


2000.79 


950 


187 


2001.16 


798 


209 


2001.70 


960 


168 



Table 4. Spiral Disk Model Parameters 



Epoch 


A 


q 


^o(°) 


1991.84 


3 


-1.0 


30 


1992.76 


3 


-1.0 


140 


1994.01 


2 


-1.0 


170 


1994.92 


2 


-1.0 


200 


1996.06 


2 


-1.0 


320 


1997.73 


3 


0.0 


420 


1998.09 


2 


1.5 


470 


1998.99 


2 


2.3 


540 


2000.79 


2 


2.7 


600 


2001.16 


2 


3.0 


640 


2001.70 


2 


3.0 


710 
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Fig. 1. — Comparison between profiles obtained at approximately the same epoch (within two 
months) but using different telescopes. Spectra are identified in the figure by the epoch (mm/yy) 
and telescope: CTIO 4m (heavy continuous line), Keck II (continuous line), KPNO 2.1m (dotted 
line), MDM for epoch 12/98 (dotted line as well), and CTIO 1.5m (dashed line). 
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Fig. 2. — Broad double-peaked Ha profiles from 1991 Nov. to 1997 Jan. (after subtraction of the 
stellar population contribution). Dotted lines show the 1991 Nov. profile superimposed on the 
subsequent profiles, for comparison. 



-23- 




Fig. 3. — Broad double-peaked Ha profiles from 1997 Sep. to 2002 Oct. (after subtraction of the 
stellar population contribution). Dotted lines show the 1997 Sep. spectrum superimposed on the 
subsequent profiles, for comparison. 
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Fig. 4. — Evolution of the measured properties of the double-peaked profiles. Left: Peak velocities 
corresponding to Gaussian curves fitted to the blue {Vb, top) and red {Vr, bottom) peaks. Right: 
Vr - Vb (top) and F^road (bottom). 
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Fig. 5. — Evolution of ratio between the blue and red peak fluxes (of the fitted Gaussians). 
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Fig. 6. — The inverse correlation between the velocity difference between the two peaks and the 
flux of the broad Ha line. 



-27- 




Fig. 7. — Fits of the elliptical disk model to the profiles from 1991 Nov. to 1997 Sep. by changing 
only (j)Q (the orientation of the disk relative to the line of sight) and (the radius of the peak of 
the emissivity law). The values of these parameters are listed in Table 3. 




Fig. 8. — Fits of the elliptical disk model to the profiles from 1998 Feb. to 2001 Sep. by changing 
only (f)o and ^q. Dotted lines show these fits blueshifted by 12 A. 
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Fig. 9. — Evolution of the elliptical disk model parameters: emissivity-law peak radius (top); 
and orientation of the disk relative to the line of sight 00 (bottom). 
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Fig. 10. — Elliptical disk cartoon, illustrating the change in the disk emissivity and orientation 
between the first and last epoch fitted with this model. Units are lO^r^ and the observer is to the 
right. Emissivity scale is arbitrary, white representing the brightest regions. 
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6400 6600 6800 

A(A) 

Fig. 11. — Fits of the spiral disk model to the profiles from 1991 Nov. to 1997 Sep. by changing (j)Q 
(the orientation of the spiral arm relative to the line of sight), q (the slope of the radial component 
of the emissivity law) , and A (the spiral arm contrast paramenter) . The values of these parameters 
are listed in Table 4. 




Fig. 12. — Fits of the spiral disk model to the profiles from 1998 Feb. to 2001 Sep. by changing 0o 
and q. The values of these parameters are listed in Table 4. 
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Fig. 13. — Evolution of the spiral disk model parameters: emissivity-law power-law index q (top), 
and orientation of the spiral arm relative to the line of sight (bottom) in the precessing scenario. 
The straight line is a linear fit to the variation of (po with time. 
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Fig. 14. — Spiral disk cartoon, illustrating the disk emissivity and orientation for the first and last 
epoch fitted with this model. Units are lO^r^ and the observer is to the right. Emissivity scale is 
arbitrary, white representing the brightest regions. 



